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We study the stability of superpositions of macroscopically distinct quantum states under decoher- 
ence. We introduce a class of quantum states with entanglement features similar to Greenberger- 
Horne-Zeilinger (GHZ) states, but with an inherent stability against noise and decoherence. We 
show that in contrast to GHZ states, these so-called concatenated GHZ states remain multipar- 
tite entangled even for macroscopic numbers of particles and can be used for quantum metrology 
in noisy environments. We also propose a scalable experimental realization of these states using 
existing ion-trap set-ups. 
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Ever since quantum mechanics was established as a the- 
ory describing light and matter at a microscopic level, the 
logical implications and its possible extension to the macro- 
scopic realm have been subject to vivid debates. Already in 
1935 Schrodinger [T] pointed out in his famous Schrodinger- 
cat gedanken experiment that quantum mechanics predicts 
the existence of superpositions of macroscopically distinct 
states, thereby highlighting the counterintuitive features of 
quantum mechanics at a macroscopic scale. While this was 
not meant as a serious experimental proposal, remarkable 
progress towards an experimental realization of such "cat- 
states" has been reported in the last decade (see e.g. [DE])- 
The interest in these states is not only triggered by funda- 
mental considerations, but also holds a more practical per- 
spective due to possible applications in quantum metrology 

SHE]- 

An archetypal "cat-state" is the multipartite 
Greenberger-Horne-Zeilinger (GHZ) state 

\GHZ±) = ±{\Q)* N ±\1)* N ). (1) 

While such a state can be argued to constitute a macro- 
scopic quantum superposition and allows for parameter es- 
timation surpassing the standard quantum limit, it is also 
widely believed that the influence of noise and decoherence 
renders the observation of related interference effects at a 
macroscopic level very difficult. It has indeed been pointed 
out with various arguments [9TTT3] that the GHZ state is 
unstable under decoherence in the sense that its quantum 
properties such as entanglement or quantum coherence as 
well as its applicability for parameter estimation vanish for 
macroscopic system size practically instantaneously, with a 
rate growing exponentially with the size of the system. 

In this letter we introduce a class of quantum superpo- 
sition states that shows similar features as the GHZ state, 
but remain multiparticle entangled even for macroscopic 
system sizes and in noisy environments. The basic idea is 
to encode the states |0) and |1) by small groups of qubits 
such that quantum properties can be conserved even under 
decoherence. This is similar to the protection of general 
quantum information using quantum error correction, how- 
ever no active correction or feedback is required. We study 
the action of decoherence on such encoded superposition 
states, which we model by single-particle white noise. We 



however point out that our observations are not restricted 
to this particular noise model, but are a generic feature of 
these states. 

We illustrate the stability of the states by calculating sev- 
eral quantities which are commonly seen as typical quan- 
tum properties and show that they are exponentially more 
stable under decoherence than for the standard GHZ state. 
These quantities include the trace norm of quantum coher- 
ence terms, entanglement measures such as the negativity 
or multipartite distillability, as well as the Fisher informa- 
tion. The latter allows us to conclude that the states in 
question can be used for parameter estimation surpassing 
the standard quantum limit even in a noisy environment. 
Finally, we propose a scalable experimental realization of 
these states using existing ion-trap set-ups. 

Definitions and notation. — Our setting is the Hilbcrt 
space of Nm two-level systems, where N is the number 
of logical qubits, each built of m physical qubits. The state 
of interest is then defined as 

\<t>c) = ^(\GHZ+f N + \GHZ-f N ) (2) 

and is called concatenated GHZ (C-GHZ), since it exhibit a 
GHZ-like superposition of a TV-fold tensor product of GHZ 
states. Note that for m = 1 it simplifies to the standard 
GHZ state up to a local basis change. We consider groups 
of size m as blocks constituting one logical qubit, and are 
interested in (entanglement) properties of the states with 
respect to these blocks. That is, we consider each block as 
one logical system on which we allow for joint operations. 

We study the stability of states under uncorrelated local 
(i.e. single qubit) decoherence processes, 

P(pc = £ l^cX^cl = £i^2 •••£jv \<t>c){<t>c\ , (3) 

where Eip = pp+^—f^ Y^j=o P a j corresponds to a single 
qubit depolarizing (white noise) channel acting on qubit i 
with error parameter p — exp(^nt), and o~j denote Pauli 
operators with o-q = 1. Notice that this uncorrelated noise 
model is rotational invariant with respect to single qubit 
operations and guarantees that no preferred basis with in- 
creased stability exists. Basis dependent noise models or 
correlated noise models typically allow to identify preferred 
basis or even decoherence free subspaces, thereby offering 
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the possibility to identify states that are stable under such 
a particular noise process. While this is certainly of in- 
terest and practical relevance in many physical situations, 
here we are not interested in such basis-dependent stabi- 
lization effects. We rather would like to identify states that 
are stable under any (local) noise process and not just for 
some special noise model |14| . 

Stability of the coherence terms. — We start by con- 
sidering the decay of coherences of the state (J3J, p G = 
£ \GHZ+){GHZ m f N . These off-diagonal elements distin- 
guish the coherent superposition of two macroscopically 
distinct states from a classical mixture, and can therefore 
be seen a characterizing the "quantumncss" of the state. 
For any state of the form \a)® N + , the trace norm 

||^4||i = tiV A^A of the corresponding coherences \a}(b\® N 
always decreases under local white noise exponentially with 
aN. However, for po we find that a can be made arbitrary 
small by increasing m. 

To calculate ||pol|ii it is sufficient to consider one 
logical qubit poi = £ \GH Z+)(GH 'Z~\, as ||po||i = 
| J pox 1 1 x • A direct calculation yields ||/9oi||i = 1 — 
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we can bound this sum from below using Stirling's formulae 
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For fixed noise parameter p, ||po||i tends to one exponen- 
tially fast with increasing m. In fact, by letting m itself 
scale logarithmically with TV, m — 0(\og(N)) we find that 
the exponential decay of ||po||i can be effectively frozen 
as limjv->oo ||po||i — ^ 1 m this case, i.e. in order to stabi- 
lize the off-diagonal elements the number of physical qubits 
per logical block has to grow only logarithmically with the 
number of logical qubits. 

Entanglement properties. — We now turn to multipar- 
tite entanglement properties of the state with respect to 
the blocks constituting the logical qubits. We investigate 
the distillability of the decohered states, as well as the de- 
cay of certain multipartite entanglement measures. A N- 
party mixed state p is called TV-party distillable entangled 
if one can generate by means of local operation and clas- 
sical communication any TV-party entangled pure state -in 
particular a GHZ state- from many copies of p |16| . The 
state p remains multipartite entangled in this sense. It has 
been shown for general noise models that the GHZ state 
looses its multipartite distillability faster with increasing TV 
|11) . leading to a lifetime of distillable entanglement that 
vanishes practically instantaneously for large system sizes. 
Similarly, for any fixed noise parameter there is a maximal 
number of blocks that can remain multiparty entangled. In 
contrast, we now show that the C-GHZ state remains TV- 
party distillable even for macroscopic iV by providing an 
explicit distillation protocol. 

First, every logical qubit of p$ c is projected into the two- 
dimensional subspace spanned by |0)® m and The 
resulting state can be viewed as a TV-qubit state consist- 
ing of the logical qubits |0 L ) = \Q)® m and \1 L ) = |l)® m . 



The state is up to a logical Hadamard operation [T7] di- 
agonal in the GHZ-basis, and its distillability can be ob- 
tained directly |16| . A measurement of all but two of the 
logical qubits in the basis {\0l) , results into a two- 

party state. The fidelity of p with respect to the (logical) 
Bell-state \<j>+) = (|0 L L ) + \1 L 1 L ))/V2 is given by F = 

<#M#) - H 1 + ^[ 1 + ( ^ 2 ] + ( ^ )Ar }' where 

d = [(l+p) m + (l-p) m ]/2 m and o = [(l+p) m -(l-p) m ]/2 m . 
As long as F > 1/2, this two qubit state is distillable, from 
which multipartite distillability follows as the results holds 
for any pair of logical qubits |16| . One observes that by 
increasing m, the range of distillability can be exponen- 
tially extended. If we consider e.g. p — 0.9, the size of 
the GHZ state is limited to TV = 53 [TT|. However, al- 
ready a block size of m — 10 allows for entanglement up to 
TV = 10 12 in the C-GHZ state, i.e. the state remains dis- 
tillable for a non-vanishing time even for this macroscopic 
system size. If we consider again blocks growing logarith- 
mically with TV, to = 0(log(TV)), we can approximate for 

small noise (f) « 1 — 2 [e(l — p)(l + ™ an d we find 
F > 1/2 for any finite TV. Notice that also for the GHZ 
state, grouping several qubits to form a block allows one 
to increase the lifetime of distillable entanglement with re- 
spect to these blocks [11], however a linear growth of block 
size, m = O(N) with respect to the number of blocks TV is 
required to stabilize the state in contrast to the logarithmic 
growth, m = O(logTV) for the C-GHZ state. 

Secondly, we consider the decay rate of entanglement, 
where we use the negativity with respect to bipartite splits 
as measure [TS] • We concentrate on the bipartition one (log- 
ical) qubit vs. rest, as this is the most fragile bipartition 
for the GHZ state [H] [H] and the C-GHZ state. Entangle- 
ment with respect to other bipartitions decays slower and 
remains non-zero for longer times. In [12 it has been shown 
that the negativity of the GHZ state shows an exponential 
decay, both with respect to time and system size N, leading 
to negligible values of entanglement for large N. Here we 
show that the decay rate with respect to N can be made 
arbitrarily small by increasing the block size m. 

The negativity of p ( f >c is defined as Af(p^ c ) = 



| nipjclll — Tr(/°^,))' wnere partial transpositions is 
meant with respect one logical qubit. By applying a logical 
Hadamard rotation |17| to each logical qubit, the resulting 
operator becomes block diagonal and can be directly diag- 
onalized. We obtain an analytic formula involving a sum of 
exponentially many eigenvalues. We find that the number 
of negative eigenvalues growth like a multinomial coefficient 
with respect to TV, m, while the modulus decreases exponen- 
tially. For system sizes up to m — 7 and N = 50, we can 
evaluate this expression numerically, see Fig. [TJi. We ob- 
serve also for large N a significant entanglement even under 
the influence of decoherence, where the negativity decays 
exponentially with N. However, the decay rate decreases 
with growing m. In fact by fitting the tail to an exponential 
function aexp(— -fN), we obtain a rate 7 that decays itself 
exponentially with m, leading to non-vanishing negativity 
even for large N. 



(a) (b) 




FIG. 1: (a) Relative negativity N(p$ c )/A/o for the bipartition 
one logical qubit vs. rest for fixed noise parameter p = 0.9 as 
a function of system size N, where the negativity of the initial 
state Ao =0.5 for all m and N. Notice the logarithmic scale, (b) 
Fisher information J~a(P()>c) f° r nxe d noise parameter p = 0.9 
and A — ^2 k (o~f m ) k , as function of system size N for different 
block sizes m. For comparison, the region between the standard 
quantum limit (SQL) and the Heisenberg limit is shaded. 

Fisher information. — We now turn to quantum metrol- 
ogy [1HE] and show that the C-GHZ states can be used for 
parameter estimation with a sensitivity above the standard 
quantum limit also in the presence of decoherence. Here a 
given quantum state p evolves under U — exp(—i9A) and is 
then measured to estimate the unknown parameter 9, where 
A is a local generating operator. The confidence interval 
for parameter estimation can be bounded from above by 
the Cramer-Rao inequality 56 > , 1 |6] , where Ta (p) 

y/ tiTa (p) 

is the so called Fisher information and n denotes the num- 
ber of repetitions of the experiment. For mixed states with 
spectral decomposition p = ^ k \k \k){k\ the Fisher infor- 
mation reads [5] 

F A {p) = 2^(A fe - A J ) 2 /(A fe + A 3 )| (k\ A \j) | 2 . (5) 

If the Fisher information is of order N -which corresponds 
to the standard quantum limit and is the upper limit for 
separable states [5]- no gain compared to classical protocols 
can be achieved. The GHZ state -as well as other entangled 
states- are however capable to reach the so-called Heisen- 
berg limit, where T A is 0(N 2 ) for the (optimal) choice 

A = Y2k=i Hi thereby leading to a quadratic improve- 
ment in the achievable precision. Under general decoher- 
ence processes the GHZ state however quickly looses its 
advantages and especially for large N there is no substan- 
tial gain as compared to uncorrelated particles |13j . 

For the C-GHZ state we show that p^, c remains for suf- 
ficiently large m above the standard quantum limit even in 
the presence of decoherence, and therefore stays useful for 
parameter estimation. We consider local generators con- 
sisting of sums of local interactions Ac — Ylk=i( a x rn ) k > 
where local is again understood in the sense of logical 
qubits. This means that for a given m, interactions involve 
blocks of m physical qubits af m . For the pure state \<f>c) 
this yields TaQ^c}) = 4/V 2 as for the GHZ state. We now 
investigate whether the decohered state p ( f >c is after some 
fixed time still useful for parameter estimation by calculat- 
ing F A {p<t> c )- After a logical Hadamard operation on each 



block [T7], the corresponding eigenvalues and eigenvectors 
can again be directly obtained due to the block-diagonal 
structure. We find that even for large N the Fisher in- 
formation stays well above the standard quantum limit, 
despite the exponentially decay with respect to the system 
size TV, see Fig. Similarly as for the negativity, the de- 
cay rate decreases with growing m. In fact by fitting to an 
exponential function aexp(— jN), we obtain again a rate 7 
that decays itself exponentially with m. 

State Preparation. — Having established that the C-GHZ 
state is robust against decoherence, we now discuss a pos- 
sible experimental realization of these states using trapped 
ions. We consider the set-up of [3]. The dominant source 
of noise in this case is correlated phase noise, which can 
be countered by using a decoherence free subspace. This 
can easily be combined with our approach by using C- 

GHZ states with logical qubits \GHZ*) = ^(|01)®" 1 / 2 ± 
|l )®m/2) for even m, thereby offering additional robust- 
ness against other noise processes as discussed above. 

Entangling operations on a linear string of ions are per- 
formed using the multipartite M0lmer-S0rensen (MS) gate 
PI, CTn(0 = I\Vi Ui= k+ i UmW with U kl (0 = e*°i*oi. 
Up to local unitaries, U n (j) transforms the state |0)® n into 
a n-qubit GHZ state. 

We now show how U n can be decomposed with single- 
qubit unitary operations [19 to obtain a C-GHZ state 
from the n-qubit GHZ state with n = mN . As a ba- 
sic tool, we use that additional single qubit operations 
on a subset of ions, Z(G) — e'^'eG 17 ! ; allow one to 
reverse some of the two-qubit interactions induced by 

u n , z({j})u n (0zt({j}) = n k <i u kl (0U k <i u kl (-0- 

k^l^j k or l—j 

In particular, applying Z(G)U n (z)Z^ (G) with G = 
{1,2..., I|} reverses all interactions between the two 
groups {1, . . . , 77} and {§ + 1, . . . , n}, but leaves the in- 
teractions invariant within the groups. An additional ap- 
plication of J7n(f ) results in an effective phase £ within a 
group, and cancelation of all phases between the groups. 
One can now apply this method inductively, by splitting 
again each group in two halves and applying Z(G') on 
G' = {1, . . . f } U {2a + 1, . . .71} while using MS-gates with 

U n (j). In total, N applications of C7jvm (iff) anc ^ N — 1 
intermediate local rotations on groups of Nm/2 ions allow 
one to obtain an operation V acting on TV blocks of size 
m, where no interactions takes place between blocks and a 
phase gate with phase | acts between each pair of parti- 
cles within each of the blocks. One can then generate first 
the GHZ state in the z-basis by applying U n (j), and from 
there the C-GHZ state by applying V. This offers a scal- 
able way of generating robust superpositions states in this 
ion-trap set-up. 

Notice that the protocol is optimal -up to at most a fac- 
tor of two- with respect the the total time the entangling 
MS-gate has to be applied, which is | in our scheme. The 
protocol is also efficient with respect to the required num- 
ber of intermediate local unitary operations and MS-gates, 
which only scales linearly with N. In fact, for small N we 
can show optimality with respect to the number of required 
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MS-gates for protocols of this kind, consisting of MS-gates 
with arbitrary £'s and local Z-gates. We remark that more 
efficient preparation schemes are possible if MS-gates acting 
on subsets of qubits are available. In particular, a scheme 
involving a total of 0{m 2 N) commuting two-qubit gates 
with interaction time f can be found. 

Discussion and Conclusion. — After illustrating the sta- 
bility of the C-GHZ state with respect to various proper- 
ties, one might wonder whether this is a generic feature of 
such encoded GHZ states, \Q L )® N +\1 L )® N . As long as |0 L ) 
and are orthogonal, the whole state keeps its "catness". 
However, we found for m — 3 and m = 4 that out of 10 6 
randomly chosen pairs of orthogonal quantum states (uni- 
formly distributed with respect to the Haar measure) , none 
showed a robustness comparable to the C-GHZ state with 
respect to decay of coherences £ |0l)(1l| 8>JV . Similarly, con- 
sidering blocks of size m for a GHZ state of size mN also 
leads to an (exponentially) smaller robustness [TT] . Perhaps 
more surprisingly, also using codewords of error correcting 
codes (a CSS code with m = 5) for |0l),|1l), we found 
that without active error correction, the coherences still de- 
cay much faster as for the C-GHZ state |2Q] . 

To summarize, we have introduced a class of quantum 



states that exhibit quantum properties similar to the GHZ 
state, but that are more stable under uncorrelated decoher- 
ence processes that we modeled by white noise. Although 
not explicitly shown here, equivalent robustness can be 
found for other decoherence models [2T]. We have demon- 
strated that the trace norm of coherences, the lifetime of 
multipartite distillablc entanglement, the decay of negativ- 
ity and the Fisher information for quasi-local generators 
can be exponentially stabilized by considering blocks of a 
small number of qubits m which are themselves prepared 
in a GHZ state. On the one hand, this may open new 
possibilities for practical applications of these states, e.g. 
in quantum communication or quantum metrology, as they 
offer an increased stability against noise and decoherence 
without need for active error correction. On the other hand, 
the existence of these states and their stability properties 
hints at the possibility to observe counter-intuitive quan- 
tum mechanical effects such as quantum superpositions of 
distinct states also at a macroscopic scale. 
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